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Abstract 

The influence of longitudinal structuring on the fast kink modes of coro- 
nal loops is investigated. Analytical dispersion relations and mode profiles are 
derived for the second-order ordinary differential equation governing the z- com- 
ponent of the perturbation in the magnetic field, SB Z . All other components 
are given in terms of SB Z . Deviations from the frequencies and mode profiles 
of homogenous loops are given as functions of the density scale height. The 
effects of the fixed and variable column masses, negative scale heights, and 
density contrasts inside and outside of the loops are studied. The frequency 
ratios, mode profiles, and the antinode shifts from those of the sine profiles of 
the homogenous loops, are tools to estimate solar photospheric parameters. To 
this end, we have expanded the relevant factors up to the second order in the 
stratification parameter. In particular, we verify that the first overtone antinode 
shifts are in the Mm range and are within the reach of the resolutions of the 
present day observations. 
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1. Introduction 

The first spatial oscillations of coronal loops were discovered in extreme ultravi- 
olet (171 A) by Aschwanden et al. (1999a). Using a MHD wave theory developed 
by Edwin & Roberts (1983), Aschwanden ct al. (1999b) and Nakariakov et al. 
(1999) interpreted them as the fundamental fast (kink) mode of coronal loops. 
Subsequently, Wang et al. (2003) observed frequency spectrum of the standing 
slow modes of the loops. 

Since then theoreticians have been supplementing the observed data, oscil- 
lation periods, loop lengths, etc., with plausible auxiliary parameters, such as 
plasma density, magnetic field strength and configuration, etc., with the aim of 
obtaining reasonably realistic models for the structure of the loops (Bennett et 
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al, 1998; Verwichte et al., 2004; Diaz et al., 2002, 2004; Van Doorsselaere et al, 
2004; Andries et al., 2005a, b, Arregui et al., 2005; Dymova & Ruderman, 2005; 
Erdelyi & Fedun, 2006; Diaz & Roberts, 2006; McEwan et al., 2006; Donnelly 
et al., 2006; Erdelyi & Verth 2007; Dymova & Ruderman, 2005; Dymova & 
Ruderman, 2006; Verth et al., 2007 Doorsselaere et al., 2007; Erdelyi & Verth, 
2007; Safari et al., 2007; Karami & Asvar 2007). 

In their recent work, Dymova & Ruderman (2005) and Safari et al. (2007) 
reduce the MHD wave equations to a single Sturm-Liouville equation for the z- 
component of the perturbation in the magnetic field. Here, we use the formalism 
of Safari et al. and show that, for an exponential plasma density stratification 
along the loop axis, the problem has a closed analytical solution. The loop model, 
equations of motion, and boundary conditions are presented in Sec. 2. The closed 
solutions, including the dispersion relation, are treated in Sec 3. Concluding 
remarks are given in Sec. 4. 

2. Description of the model and equations of motion 

A coronal loop is approximated by a cylinder of length L and radius R. Loop ends 
are fixed at the photosphere. Loop curvature is neglected, on account of R « L. 
No initial flow is assumed inside the loop. A uniform magnetic field along the 
axis pervades the loop, B = Bz. Gas pressure, gravity and all dissipative and 
viscous forces arc neglected. The density is discontinuous on the lateral surface 
of the cylinder and varies exponentially along the axis. Thus, 

p(r,z,e) = pi(e) exp(— ez/L), < z < L/2, inside tube, 

= p e {s) exp(— sz/L), outside tube, (1) 

where e is the density scale height parameter, and pi, p e are the interior and 
exterior footpoint densities, respectively. The assumption of exponential density 
is in accord with the findings of Aschwandcn et al. (1999). They conclude this 
from their stereoscopic analysis of 30 loop oscillations in EUV. Restriction of z 
to the interval [0 — L/2] is permissible on account of the symmetry of the loop 
configuration about its midpoint 

The linearized ideal MHD equations are 



| + V.(Mv) = 


0. 


Continuity Equation, 


(2) 


^ = — (Vx 
dt 47r/ 


<5B) x B, 


Momentum Equation, 


(3) 


dSB „ 

dt =VX( ' V 
V.5B = 0, 


xB), 


Induction Equation 
Solenoidal Constraint, 


(4) 
(5) 



where <5v and (5B are the Eulerian perturbations in the velocity and magnetic 
fields, respectively. An exponential <p and t dependence, is assumed, exp[— i(mcp— 
cut)] . By straightforward calculations one can express all components of Sv and 
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in terms of SB Z . The latter, in turn, is obtained from the following second 
order PDE (See Safari et al. 2007), 

2 Id d 2 m 2 a, 2 \ SB Z 
or z r or oz 2 - r z v A ) B 

where va{z) = B / y/Anp(z), local Alfven speed, has different values inside and 
outside of the loop. Equation (6), admits of a separable solution 8B z /Bq = 
R(r)Z(z), where R(r) satisfies Bessel's equation and will not be further referred 
to here, and Z satisfies the following 



d 2 Z(x) 



fl 2 e- ex Z(x) = 0, < x = z/L < 1/2, (7) 



n 



dx 2 

2 L 2 uj 2 ft(e)+p e (e) 



where 51 and pi_ e — Pi, e (e) / Pi.e(0) arc dimcnsionlcss frequency and footpoint 
densities, respectively. 

Equation (7) is an eigenvalue problem weighted by exp (—ex). Changing the 
variable xtofl cxp(— ex/2) reduces Eq.(7) to a Bessel equation with the following 
solutions 

Z(x) = ci J (2-e--/ 2 ) + C2 y (2-e— / 2 ), (8) 
e e 

where c\ and ci are constants, and J and Y are Bessel functions of first and 
second kind, respectively. The boundary conditions are 

for odd modes j ^(iJ^L ( 9 ) 
for even modes < ^z(l/2) — 
Imposing the boundary condition at z = on Eq. (8) gives 

Cl =Y (2-), c 2 = -Jo( — ) (11) 

e e 

Substituting these coefficients in Eq. (8) and imposing the boundary conditions 
at x = 1/2 gives the dispersion relations 

-r (2-)Ji(2-e- £/4 ) + J (2-)Y 1 (2-e- £ / 4 ) = 0, odd modes, 

£ £ £ £ 



-Y (2-)J (2-e- £/4 ) + Jo(2-)ro(2-e- e/4 ) = 0, even modes. 

£ £ £ £ 



(12) 
(13) 
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Figure 1. Fundamental and overtone frequencies versus e. Solid lines for fixed column mass 
but variable footpoint densities, dashed lines for variable column masses but fixed footpoint 
density. All frequencies are in units of ttua (e = 0)/L. 




Figure 2. Mode profiles, Z n (z), a ,b , c, d, corresponding to n = 1,2,3,4, respectively. Solid, 
dashed, and doted lines are for e = 0, 2, and 10, respectively. 



Equations (12) and (13) are similar to those of Diaz & Roberts (2006) with W — ► 
in their analysis. In the remainder of this section they are solved analytically 
for weakly stratified loops and numerically for arbitrary stratifications. 
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2.0.1. Weak stratification 

Bessel functions can be expanded as 



J v (z) = \ — {P„(z) cos z' - Q v (z) sin z'} , z' = z - ( v + \ J ^, 

V 1TZ V I J I 

Y A Z ) = \ —{Pv(z) sin z' + Q v (z) cos z'}, 

V TTZ 



(14) 



where 



(4* a -l)(4^-9) 1 (4^»-l) 1 

128z2 + L V ; ' 8z + V 



We insert Eq. (14) into Eqs. (12) and (13) and reduce them for e <C 1 and 
find 



odd, 



2n(l + e^+£ 2 (|.- even, 

(15) 

R — - 

where n = 1, 2, • • •, and wi = [27r(/jj(0) + p e (0))] 2 is the fundamental kink 
frequency of homogenous loops. Expectedly, w n — > noji as e — > 0. The following 
ratios are noteworthy: 

lj 2 _ i _ „ 101 i P 2 1 2525 , 3535 \ 

2o)i ^TStp 3 ' ' fc V64tF " r 2048tt 2 I 

25 (2n-l) 2 -l _ 2 125 I 5,7 \ (2n-l) 2 -l 



(2n-l)wi 



— 1 e- 25 l-m-i; -l , ,-2 125 I 5,7 \ \tn-V) -l /-. fi x 

- 1_e 4^ (2n-l) :i +£ W^\^ + W> (2n-l)'> l ib J 



c^2n it ,-2 1 n — 1 

— -L ~T~ c — T 3 — 

nuJ2 lb7T z 

The frequencies and the ratio of any two odd numbered frequencies begin de- 
creasing linearly with e. The ratio of two even modes, however, begins decreasing 
quadratically with e. These features are also seen on the diagrams of Fig 36. 
Observational verification of these points, however, has to await the availability 
of more extended and higher resolutions data. Presently only two frequencies in 
three loops are available ( Verwitche et al. 2004, Van Doorsselaere 2007). 

2.0.2. Arbitrary stratification - Numerical approach 

We use Newton-Raphson's numerical method to solve Eqs. (12) and (13) for the 
eigenfrequencies. In the range, < e = L/H < 20, the fundamental and three 
higher kink frequencies, w„, n = 2, 3, and 4, and the ratios, u) n /u>i, are computed. 
The data are plotted in Fig. 3 for two density contrasts, p e (e)/pi(e) = 0.1 and 
0.5. As the density contrast increases the frequencies shift down. Their ratios, 
however, remain unchanged. The ratio u] n /u)\ begins with n and decreases with 
increasing e, in compliance with Eqs. (16). 
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Andries et al. (2005b) maintain that the frequency ratios could be used as 
a seismological tool to estimate the coronal density scale heights. From the 
TRACE data, Verwichte et al. (2004) find the ratio u 2 /ui to be 1.64 and 1.81 
for loops designated by D and C in their records, respectively. Van Doorsselaere, 
Nakariakov, & Verwichte (2007) revisited the same ratios from the observational 
data to be 1.58 and 1.82 for the same loops, respectively, and 1.795 for another 
loop in their current analysis. With the help of Fig. 1, we find the corresponding 
e to be 4.98 and 1.9, respectively. Assuming typical loop lengths, L = 100 — 400 
Mm, the density scale heights fall in the range of H = e^ 1 L ~ 20—82 and 53—210 
Mm, respectively. These scale heights are slightly different from the findings of 
Andries et al. (2005a, b) and Safari et al. (2007) for a sinusoidal density profile. 

A noteworthy point is the effect of column mass on frequencies. In Fig. 1 
we assume a constant footpoint density contrast, p e / Pi = 0.1, and vary e. Con- 
sequently, the total column mass of the loop changes. Compared with variable 
density contrast but fixed column mass, the mode profiles and the frequency 
ratios remain unchanged. The frequencies themselves, however, behave slightly 
differently. For variable column mass models, the frequencies increase more 
sharply with e. 

The mode profiles, Z{z) of Eq. (8), are shown in Fig. 2 for n = 1,2, 3, and 4. 
For the unstratified case, e = 0, the profiles arc sinusoidal. With increasing e, 
they depart from the sine curves. Antinodes shift toward the footpoints. Stronger 
the stratification, the greater the shift is, in agreement with the findings of Safari 
et al. (2007) and Verth et al. (2007). Van Doorsselaere et al. (2007) point out, the 
shift of the antinodes is potentially, a coronal seismological tool to estimate the 
density scale heights. In Fig. 3, we have plotted the antinodc shift, zf n — zf™ , of 
the first overtone versus e. The shift w 0.02e, grows approximately linearly with 
e. Our numerical result shows that, the shift in the antinode for different density 
contrasts , p e /pi = 0.1 and 0.5, are the same. For typical loops, of lengths 
100 — 400 Mm and density scale heights, H = 50 and 100 Mm, the antinode 
shift falls in the range 2.85 < z^ n - zfj? < 56.64 Mm and 1.35 < z^ n - z^ < 
25.12 Mm, respectively. Observation wise, the resolution of current solar satellite 
facilities, e.g., TRACE, SDO, SO, etc., seems adequate to detect such antinode 
shifts and estimate the density scale height of solar coronae. Verth et al. (2007) 
study semi circular loops of sinusoidal density profiles and find the antinode 
shifts 0.028e. 

Our numerical results show that, for a given e, the mode profiles arc insensitive 
to changes in density contrast, p e (e)/pi(e). The differences between the mode 
profiles of the stratified and unstratified cases, AZ n = Z n (e,z) — Z n (e = 0, z), 
n = 1,2,3, and 4 are plotted in Fig. 4. Erdelyi & Verth (2007) maintains that 
these differences in the mode profile are so small to be resolved by current EUV 
instruments of TRACE mission. 

Another interesting point is negative scale heights, suggested by Andries 
(2005a, b) on the basis of the error bars in the observations of Verwichte et 
al. (2004). Here, the density at the apex is higher than at footpoints. Unlike the 
positive scale scale height case: a) u>\ decreases with increasing |e|, see Fig. 5; b) 
higher overtones, uj n , n = 2, 3, •• •, however, increase with e|, though at slower 



anal_sol.tex; 10/02/2008; 14:08; p. 6 



Exact solutions of standing kink modes 




e 



Figure 3. Antinode shift, z^" — z^™ , (normalized to L) against e. The shift varies almost 
linearly with e 

(a) (b) 




Figure 4. Differences between the eigenprofiles of the stratified and unstratificd cases, AZ[, 
are plotted z. 



rate; c) the ratios u) n /w\ increase with |e| (not presented in a diagram). As |e| 
increases, the mode profiles and their node and antinodes move away from the 
footpoints and concentrate more and more in the inner regions of the loop. 

3. Conclusions 

Suggested theoretical models of 3D coronal loops are, still, far from the realities. 
Many complicating factors, such as variable cross sections, variable magnetic 
fields, non-zero (3 plasmas, etc., are to be accounted for in a realistic study of 
both the equilibrium structure and the perturbed state of actual loops. Here, 
we study the oscillations of loops with exponential density variations along the 
loop axis. 
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Figure 5. For the case of negative scale heights, the frequencies plotted versus e, for two 
different density contrast, pe/pi = 0.1 (solid lines) and pe/pi = 0.5 (dotted lines). 



- Analytical dispersion relations, Eqs. (12) and (13), and analytical mode pro- 

files, Eq. (8) and Fig. 2, are derived. 

- For weak stratifications, the kink frequencies and the frequency ratios are found 

up to the second order in e, Eqs. (15) and (16). 

- Increasing the density contrast decreases the frequencies but their ratios and 

shape of the profiles remain unchanged. 

- Models with variable total column mass, but constant footpoint densities, are 

investigated. Compared with models of constant total mass, the frequencies 
increase more sharply with increasing e. 

- The case of negative scale heights is investigated and results are compared 

with those of positive e's. 

- For 1.58 < ui/w-i < 1.82, and for typical loop lengths, 100-400 Mm, the density 

scale heights fall in the range of 20-210 Mm, in agreement with Andries et 
al. (2005a, b), Safari et al. (2007), McEwan et al. (2006), and Donnelly et 
al. (2006). 

- Based on our simple theoretical model and typical coronal conditions, the 

antinodc shift of the first overtone mode profiles are in the range of 1.3—56.6 
Mm. They arc in the range of the detectability of the resolution of the 
current observational instruments. 
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